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GRAVITY AIDED INERTIAL NAVIGATION 
SYSTEM 


FIELD OF THE INVENTION 


This invention relates to inertial navigational systems. 
More specifically, it relates to a new integration of a 
conventional three axis Inertial Navigation System 
(INS) with gravity gradiometer instrument(s) and a 
height or depth sensor. 


BACKGROUND OF THE INVENTION 


Inertial navigation systems employ a triad of acceler- 
ometers in combination with a triad of angular rate 
sensors for providing velocity and position information 
to a navigational computer. Such systems provide con- 
tinuous measures of vertical acceleration, (perpendicu- 
lar to the earth's surface) and horizontal acceleration, 
(parallel to the earth's surface) and can determine verti- 
cal and horizontal velocity as well as vertical position 
and horizontal position by time integration of the verti- 
cal and horizontal acceleration. 

Simple integration of vertical and horizontal acceler- 
ation to provide vertical and horizontal velocity, and 
integration of the velocities to provide vertical and 
horizontal position results in unbounded errors in the 
resultant vertical and horizontal axis information due to 
two causes. First, INS accelerometers measure the de- 
sired vehicle acceleration but in combination with the 
gravity field. Although the normal component of the 
gravity field can be accounted for, the anomalous com- 
ponent cannot. Consequently, this acceleration together 
with accelerometer and gyro instrument errors lead to 
navigation errors at Schuler and siderial frequencies 
which grow unbounded with time. Second, any vertical 
position error results in an error in the computer de- 
rived normal gravity component which causes the ver- 
tical position error in the INS to increase without 
bound. 

A gravity gradiometer measures the gravity field 
independently of vehicle acceleration so it can be used 
to correct INS accelerometer outputs for the anomalous 
gravity field component. The unstable INS vertical 
channel is often remedied in practice by integrating a 
height sensor (depth gauge, altimeter, or surface ship 
sea level knowledge). The difference between inertial 
vertical position and that of the height sensor is used to 
stabilize the inherently unstable vertical channel. 

The resulting navigation performance, although 
much improved over that of the conventional unaided 
INS, continues to exhibit growing errors at schuler and 
siderial frequencies due to INS instrument (including 
gravity gradiometer) errors. 

If a gravity field map is available the integrated INS 
(conventional INS-r-gravity gradiometer sensor(s)- 
+ depth sensor) can implement a map matching mode. 
In this scheme mapped values of gravity field anomaly 
are compared with measured values in a filter to bound 
velocity and position errors. 

Apart from the requirement for a gravity map, navi- 
gation performance is somewhat compromised because 
the core INS around which map matching navigation is 
implemented has unbounded velocity and position er- 
rors so in weak gravity signature regions navigation 
errors will tend to grow. 

Although the prior art which integrates a conven- 
tional 3 axis INS with gravity gradiometer(s), and a 
height sensor greatly improves navigation performance 


10 


15 


25 


30 


35 


40 


45 


50 


55 


60 


65 


2 


in both the INS (no gravity map) and in the map match- 
ing mode, it does not fully utilize overall navigations 
system velocity error observability. If this velocity 
error observability is exploited these errors can be 
bounded even in the INS mode. The resulting improve- 
ment in the core INS mode, where now only east posi- 
tion error is unbounded, leads to improved map match- 
ing performance. 


OBJECTIVES OF THE INVENTION 


A primary objective of the invention is to provide an 
autonomous covert Inertial Navigation System (INS) 
wherein inertial velocity errors are bounded without 
external navigation aids or active instrumentation of 
ground speed. 

An objective of the invention is to provide accurate, 
self-contained, passive navigation with bounded veloc- 
ity error and improved position error (bounded Schuler 
and siderial components) by integrating gradiometer 
and height (depth) sensors with conventional inertial 
navigation instruments. 

A further objective of this invention is an integration 
Scheme which takes advantage of navigation system 
velocity error observability wherein east velocity error 
through the mechanism of vertical Coriolis accelera- 
tion, is manifested in observable height (depth) error. 

A still further objective of this invention is an integra- 
tion scheme wherein gravity disturbance vector esti- 
mates, based mainly on gradiometer measurements, are 
used to compensate inertial navigation system acceler- 
ometer measurements. 

A further objective of this invention is an integration 
scheme which takes advantage of navigation system 
velocity error observability wherein erroneous velocity 
estimates used to propagate gravity states cause gravity 
State errors which are manifested in observable height 
error and gradiometer innovation errors. 

A still further objective of this invention is an integra- 
tion scheme resulting in a core INS system with all 
velocity and position errors bounded except east posi- 
tion which may be controlled by a star tracking or by a 
map matching mode. 


SUMMARY 


Integration of gravity gradiometer sensor(s) and 
height (depth) sensor with a conventional inertial navi- 
gation system (INS) results in a new INS termed the 
Gravity Aided Inertial Navigation System (GAINS). 
Like the conventional INS, GAINS is autonomous and 
covert but additionally has bounded velocity error re- 
sulting in bounded Schuler and sidereal errors. Gradi- 
ometer(s) and height sensor measurements differenced 
with predicted values are processed by the system to 
obtain estimates of the anomalous gravity field, naviga- 
tion errors, and instrument errors. 

One open loop and three closed loop mechanisms are 
used to control velocity error. First, real time estima- 
tion of the gravity field includes disturbance vector 
estimates which are used to open loop compensate the 
INS accelerometers so that navigation errors are not 
driven by the disturbance vector. Second, subtracting 
gravity field estimates from accelerometer measure- 
ments results in estimated vehicle accelerations consist- 
ing of true vehicle acceleration components and veloc- 
ity dependent Coriolis terms. After using INS estimates 
of velocity to compensate for Coriolis effects, accelera- 
tion estimates are integrated into vertical and horizontal 


5,339,684 


3 


estimates of velocity and position. The vertical position 
estimates are differenced with height sensor measure- 
ments to obtain an observation of vertical channel esti- 
mation errors. Since East velocity error leads to Corio- 
lis compensation error in the vertical channel, it be- 
comes observable by the optimal filter. This is one of 
three closed loop mechanisms for controlling velocity 
error. 

The second and third closed loop mechanisms for 
velocity control are as follows. The gravity field esti- 
mates are propagated forward at seconds using the 
current best estimate of velocity. The linearized rela- 
tionship of gravity disturbance error to velocity error is 
given by the product of At and the estimated gradients. 
Since the vertical component of gravity disturbance is 
used to compensate the vertical accelerometer, this 
velocity dependent error is observable in the height 
error signal. The difference between predicted and 
measured gradients are likewise attributable to velocity 
and gravity field estimation error. In this case the linear- 
ized relationship to velocity error is the product of at 
with the third order gravity field estimates. 

Thus GAINS, in the pure INS mode (no gravity 
map) provides bounded velocity error and hence, 
bounded schuler and siderial errors. In this mode only 
east position error is unbounded. In order to bound east 
position error as well, a star tracker can be used or if a 
gravity map is available, a map matching outer loop can 
be implemented. 
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The Gravity Aided Inertial Navigation System is: E 


self-contained (no external signals required); covert (no 
signal emanated; near surface operation not required).; 
unrestricted in operating area; able to operate in quiet or 
active gravity regions; and able to be configured with 
proven instruments. J 
The instrument which enables this type of operation 
is the gravity gradiometer instrument, developed by 
Bell Aerospace Textron for the Unisys Corporation and 
the Strategic Systems Program Office (SSPO). This 
instrument and gravity meter assemblies were incorpo- 
rated in the Trident 11 Gravity Sensors System (GSS). 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a conceptual drawing of the rotating gravity 
gradiometer instrument. 

FIG. 2 is a cut-away drawing illustrating the physical 
orientation of the major elements of the gravity gradi- 
ometer instrument (GGI). 

FIG. 3 is a schematic illustrating the umbrella orien- 
tation of three GGIs where each spin axis has the same 
orientation with respect to the local gravity vector. 

FIG. 4 is a cut-away drawing illustrating the position- 
ing of the GGIs within an inertial platform. 

FIG. 5 is a block diagram of a basic Gravity Aided 
Inertial Navigation System. 

FIG. 6 is a block diagram of a basic Gravity Aided 
Inertial Navigation System augmented with map 
matching. 


DETAILED DESCRIPTION OF INVENTION 


To clearly present the concept on which the inven- 
tion is based, gravity field nomenclature and definitions 
are standardized in the following specification. In all the 
descriptions that follow a carat ( ^ ) is used to denote the 
best estimate of some system parameter (i.e. Vy means 
best estimate of North Velocity). A tilde (—) is used to 
denote the error in that estimate. Throughout the speci- 
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fication and claims, a height sensor is synonymous with 
depth sensor and both are vertical displacement sensors. 


GRAVITY FIELD NOMENCLATURE AND 
DEFINITIONS 


Gravitation is a potential field defined by a scalar 
potential wherein the anomalous component of the field 
is given by; 


T=W-—U 


where: 
W —total gravitation potential 
U=regular potential associated with an idealized 
regular earth 
Spacially differentiating the anomalous potential with 
respect to north, east, down (NED) coordinates results 
in the anomalous disturbance vector designated T; 


aT 


a 
T; = (TNTETD) = = 3r. 


aE 


' denoted transpose. 
A single axis accelerometer measures one component of 
the specific acceleration given by; 


Differentiating the disturbance vector spatially re- 
sults in the anomalous gravity gradient second order 
tensor designated Tj 


TNN TNE TND 


Tj = | TEN TEE TED |= 
TDN TpE TDD 
aT T T 
ane NoE ND 
2T BT aT 
EN E 2232 
2T 2T o $T 
DaN DE ap 


This tensor is symmetric since the order of differentia- 
tion does not affect the result. Since Laplace's con- 
straint applies, the tensor has five independent elements. 
A gradiometer measures components of Wj— (Uj4- T j), 

Continued spatial differentiation is used to compute 
the higher order anomalous gravity tensor terms if they 
are required. The third order term Tj for example, is 
used in gravity gradiometer map matching. 


ROTATING GRAVITY GRADIOMETER 


FIG. 1 depicts as by way of example the Bell Aero- 
space rotating gravity gradiometer instrument (GGI) 
which consists of four single axis accelerometers (A1, 
A2, A3 and A4) mounted tangentially at cardinal points 
on a rotating instrument block or wheel. Two outputs 
of the GGI are obtained by demodulating the signal 


AG=A+A2—A3—Ag 
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at sin 20t and cos 20t where Q is the wheel rotation 
rate (x} Hz). Note that if the accelerometers are per- 
fectly aligned, scale factor balanced, and linear, no 
wheel acceleration or angular acceleration about its 
spin axis appear in Ag. In addition any residual linear 5 
acceleration sensitivity will be modulated at Q and, 
therefore, will not appear in the 20 demodulation out- 
puts. So the gradiometer which is mounted on a stabi- 
lized platform is insensitive to vehicle acceleration to 
first order. The instrument measures the components of 10 
the gradient lying in the plane of the wheel. The gravity 
field at any point in be plane of the wheel is given by the 
Taylor series expansion about the wheel center. 


7)۵ = (jo + (Todi + lr Urgdodid; ^... 5 
The output of accelerometer A is 
A1-[40— TRY 1 20 
where Rj— radius vector to A1 
where A,—acceleration of the wheel center. 
and 11, is a unit vector defining the sensitive axis of A. 
To first order 25 


TxxTxy || موم‎ —sinQt 
Ab [^ 70 - Pq `L eost 


Since Ac=Aír+A2—As—A4, the first two terms in 30 
(11) will cancel while the third term will add yielding. 


š 35 
cai TxxcosQt TyysinQt .[ -singe ۶ 
Tosh! TysinQr cosQt 
Txx — T, 
AG — 4R |= > sin2Qr + 1 7 


Sin 20 and cos 20t demodulation of Ag and low pass 
filtering produces the inline and cross outputs 
2R(Txx—Txy) and 4RT respectively. 

FIG. 2 illustrates the major components which com- 45 
prise the gravity gradiometer instrument or GGI. The 
system includes four accelerometers positioned at 90° 
increments tangentially around a rotatable accelerome- 
ter support mounting block 33. A shake actuator 31 and 
shake PO accelerometer 32 are positioned along the 
spin axis of the mounting block 33 which is supported 
by axial shake flexures 34 and an axial shake adjustment 
spring 35. 

Each of the four accelerometers 21 is contained 
within its own magnetic shield housing 23 along with its 
accelerometer preamplifier 22. Additional electronics 
are located on printed circuit boards 24 which are 
mounted to the accelerometer mounting block 33. 

The accelerometers 21 and 32, shake actuator 31 and 
printed circuit boards 24 are contained within a mag- 
netic shielding housing 25 which is supported for rota- 
tion about the spin axis along with the accelerometers 
and other components mounted on the rotatable accel- 
erometer mounting block contained therein. The rotat- 
ing housing is supported by ball bearings 36 and driven 
by a motor 37 which spins the rotatable housing and the 
rotatable support mounting block contained therein 
within the relatively stationary housing 41 which in- 
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cludes a further magnetic shield 42. Slip rings 38 pro- 
vide electrical connections required to drive the rotat- 
ing elements and power the subsystems and extract the 
data. Rotational speed is controlled via a shaft encoded 
disk 26 and encoder pickoff 27. 

The gravity gradiometer instrument component 
housing 41 includes a mounting plate 43. The mounting 
plate is used to firmly secure the gravity gradiometer 
instrument within the inertial platform as illustrated in 
FIG. 4 wherein two of the three gravity gradiometer 
instruments 44 are depicted. 

Since the gradient tensor has five independent com- 
ponents and a single GGI measures two linear combina- 
tions in the wheel plane, it follows that if the complete 
tensor is to be measured, three GGIs are required. 

The three GGIs are relatively arranged in an um- 
brella orientation as illustrated in FIG. 3 where each 
spin axis has the same orientation with respect to the 
local gravity vector. Two of the GGIs 44 are visible 
through the cutaway illustration of the platform in FIG. 
4. Note that this cutaway also indicates the pitch, azi- 
muth and roll axis of the platform relative to the GGIs 
and the positioning of the gyros 51, horizontal platform 
accelerometers 52 and gradiometer (vertical accelerom- 
eter) 53. 

It is not necessary in all applications to instrument the 
complete tensor. For instance, a parametric study may 
be accomplished for the performance of a single, double 
or triple GGI configuration. In any event, the entire 
stabilized platform which carries the INS sensors as 
well as the gravity sensors is carouselled about local 
vertical at 500°/hr. Consequently, a vehicle moving 
slowly through a slowly changing field with a single 
gradiometer at the umbrella angle will in time measure 
the full gradient. i 


GAINS SYSTEM CONCEPT 


FIG. 5 is a block diagram of a GAINS implementa- 
tion. It has been greatly simplified to clarify the essen- 
tial features of the concept. The system consists of a 
gradiometer, a three axis inertial navigation system 
(INS), a height sensor, and an optimal filter whose func- 
tion is to integrate the subsystems to produce the best 
possible navigation (i.e. provide the best estimates of 
INS position and velocity errors, gravity field, and 
instrumentation parameters). Simply put, the invention 
is based on the observation that over time, almost all 
system errors that impact inertial navigation accuracy 
manifest themselves in the vertical channel of the navi- 
gator. If an independent measure of vertical navigation 
is available, it can be compared with the INS vertical 
channel and the difference processed to provide better 
estimates of the system errors. Different INS configura- 
tions (e.g. local-level, space stabilized, strap down) lead 
to different characterizations of system dynamics. A 
model of these dynamics is required in applying this 
invention to a specific configuration. Also open to 
choice is the type of filter used to process measurements 
(e.g. Kalman fitter, non-linear estimator, least squares 
estimator). The crucial element is knowledge of how 
overall system errors manifest themselves in the mea- 
surements. 

FIG. 5 shows the vertical accelerometer of the INS 
separately to emphasize how instrument corrections are 
dependent on estimated position and velocity. 

Included in the optimal filter is a mechanization of the 
dynamics of the INS including standard siderial loop, 
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Schuler and vertical channel mechanisms. In this way 
the propagation of system errors in the navigation chan- 
nels can be tracked. Of special interest is the depen- 
dence of both Coriolis compensation error and vertical 
gravity estimation error on velocity error. These errors 
propagate into vertical position error and thus a com- 
parison of INS vertical position with height, sensor 
measurements allows estimation of these INS velocity 
errors and all other system errors that lead to velocity 
errors. 

The dependence of gravity estimation error on navi- 
gator velocity error is implicit in any mechanization of 
gravity estimation that includes gradiometers. The fol- 
lowing is an example of this dependence using a specific 
statistical representation of the anomalous gravity field 
and its modeling as a markov process appropriate for 
state space implementation. 


A STATISTICAL MODEL OF THE 
ANOMALOUS GRAVITY FIELD 


A number of different statistical characterizations of 
the gravity field have been developed. The Sperry 
Three Dimensional Analytic Gravity Model (STAG) is 
presented as exemplary. The STAG gravity potential 
autcorrelation function is defined as: 


Ten‏ ھ2 
N= 2 CD = Ta oa‏ 


Each of the n terms covers a different range of spatial 
frequencies. Empirical fits (o-7, an, dn) to ocean survey 
data have been computed using typically N —3 or more 
terms. Statistical characterization of higher order tensor 


=. 200) 
(Grlx = [>° we 


terms and cross correlations are computed by spatially 
differentiating this equation. 

The cross correlations between the N gravity pro- 
cesses are assumed to be zero. The gravity field at a 
point can then be characterized by the sum of N gravity 
processes Gy where 


EXAMPLE MARKOV MODELING OF THE 
GRAVITY FIELD 


The STAG model of the gravity field is not Markov 
and, therefore, not suitable for incorporation into a state 
space Kalman filter. If implementation of such a filter is 
desired, a recursive model with statistics approximating 
those of the STAG model is required. In particular, a 
transition matrix , and a plant noise matrix Q, are 
needed to characterize the process noise for each G; 
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process. For a given change in position, Y, we can 
approximate the optimal transition matrix , by 

«(Y = C40.4654; + Y) C, 0.46547) Eq. 1 

اجس =- 

where d, has the same direction as Y and magnitude 
equal to the STAG parameter dy. In order for the statis- 
tics of the process to remain stationary the Q, matrix 
must have the form 


Qn=Cul0)— 643) 0070 (Y) 


This Q, is associated with plant white noise driving only 
the highest order tensor terms in Gy. This formulation 
has been found to give a markov gravity model having 
good statistical agreement with the STAG model. 


EFFECT OF VELOCITY ERROR ON THE 
FILTER TRANSITION MATRIX 


The matrix ہر‎ determined above is correct only if 
the true change in position, Y, is precisely known 
(Y «V At where At=iteration interval). In fact, this 
vector will be in error by VAt. Hence, the proper transi- 
tion matrix must include the propagation of velocity 
errors into gravity errors. In the transition matrix given 
in Eq. 1, gravity estimation errors would appear to be 
the result of previous gravity errors and driving noise 
only. In a transition matrix that accurately reflects the 
effects of velocity error, three more columns of Pn, 
need to be computed, namely those multiplying the 
three navigation velocity error states. In essence these 
error propagation terms come from integration of the 
gradients over an incorrect vector position change. In 
general 


Equation (2) 
a اک لٹ‎ 
a 36 ۲۸۷۵٢ 
(Gn)x—1At PA (Gn)x—14t AY (Gi) k-1At ای‎ 
Vpht 


For the given markov model, Eq. 2 is the required 
formulation of gravity estimation transition reflecting 
the dependence upon velocity error. 


EXAMPLE FILTER IMPLEMENTATION 
(KALMAN FILTER) 


One possible implementation of this invention uses a 
Kalman filter as the parameter estimation routine. Com- 
bining the state space modeling of the gravity field 
transition given above with standard INS state space 
models gives the system formulation required by the 
Kalman fitter. In this implementation, INS accelerome- 
ters are compensated for Coriolis effects (using naviga- 
tor velocity estimates), normal gravity (using navigator 
position estimates), anomalous gravity (using filter 
gravity estimates), and instrumentation errors (using 
filter estimates), and these accelerations are processed 
to give inertial velocity and position estimates. The state 
vector (estimates of gravity field and system parame- 
ters) and state uncertainty matrix are transitioned in the 
standard fashion to the new time and place as estimated 
by the INS navigator. A measurement vector is then 
formed. The first measurement is the difference be- 
tween the INS estimate of vertical position and the 
estimate from the height sensor. The other measure- 
ments are the gradiometer outputs. It is important to 
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note that the gradiometer need not measure the com- 
plete gravity tensor for the GAINS concept to be effec- 
tive. 

The standard Kalman gain matrix is then calculated. 
The gains are applied to the difference between the 
measurement vector and the predicted measurements. 
The results are used to improve the system estimates, 
and the state uncertainty matrix is updated as is standard 
in a Kalman filter. This procedure is repeated at a rate 
mostly dictated by the speed of the computer hardware 
chosen for a given system. 

When this system reaches steady state, all navigator 
errors will be bounded except for east position errors of 
near zero frequency (non-sidereal, non Schuler errors 
associated with polar axis drift in inertial coordinates) 

Thus a method of inertial navigation is effected 
where the velocity error and position error are bounded 
at Schuler and sidereal frequencies without external 
navigation aids or active instrumentation of ground 
speed by optimally integrating available measurements. 
For instance, bounding the velocity error and position 
error may be accomplished by comparing gradiometer 
measurements with predicted gradiometer outputs and 
combining this with knowledge about the dependence 
of this difference on total system errors. Additionally, 
this bounding of the velocity and position errors can be 
accomplished by comparing an estimated depth ob- 
tained from inertial instruments with an estimate from a 
height (depth) sensor and combining this with knowl- 
edge about the dependence of this difference on total 
system errors. 

A navigation system with all errors bounded can be 
configured using the GAINS concept as a base and 
augmented by additional information or measurements. 
An example of GAINS augmentation that does not 
require any additional instrumentation is the inclusion 
of gravity map matching into the basic system for 
bounding east error by gravity field map matching. 


EXAMPLE AUGMENTATION, GRAVITY MAP 
MATCHING WITH GAINS 


GAINS used with no maps, controls navigation er- 
rors very well with only the east position error being 
unbounded. Consequently the system has excellent 
gravity map matching performance. FIG. 6 illustrates 
the GAINS map matching block diagram. 

In the following discussion, it is assumed that opera- 
tionally the map has been preprocessed to provide grid- 
ded consistent values of T, T; and Ty at mean sea level 
with algorithm capability to interpolate and down- 
ward/upward continue the map values to the GAINS 
system location and to also provide (Tj) map data. 

In the map matching mode, the GAINS optimal filter 
carries a Markov gravity model of the map error and 
not of the full anomalous gravity field. If for example 
the map used is based an GEOSAT data then the map 
error model must characterize the higher spatial fre- 
quencies not present in the map. The INS accelerometer 
and gravimeter measurements are now compensated 
with both map values and map error estimates of distur- 
bance vector components. The gradiometer measure- 
ments are likewise compensated with both map and map 
error gradients. 

The optimal filter requires a linearization of naviga- 
tion errors appearing in the gradiometer innovation. In 
the NED frame the associated gradient error is given 
by: 
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This gradient (a function of position error P) is resolved 
into the instrument frame of each gradiometer and the 
associated in-line and cross gradient contributions are 
determined. 

The GAINS concept combines gradiometer and 
height sensor data with a three axis INS to produce a 
totally autonomous and covert navigation system hav- 
ing almost all errors (East position error excepted) ac- 
curately bounded. Inclusion of an a priori map bounds 
all errors. Since a GEOSAT map is available for all the 
oceans, the system can be run continuously in the map 
matching mode. 

Mission scenarios requiring precise relative naviga- 
tion (as opposed to absolute navigation) are satisfied by 
the GAINS system. An example of such a mission 
would have an AUV enter a poorly mapped area, 
search for and map mine locations then return to and 
rendezvous with the deploying submarine. The AUV 
simultaneously maps the gravity field of the area while 
establishing mine locations in this map frame. Any other 
vessel equipped with this map and GAINS can enter the 
area and position itself accurately with respect to the 
mapped mine field. 

Such mission scenarios and system configurations and 
operations are presented herein are by way of explana- 
tion and in no way to be considered limiting. Further- 
more many other embodiments which will readily be 
apparent to those skilled in the art may be made without 
departing materially from the spirit and scope of this 
invention. 

What is claimed is: 

1. A navigation system, comprising: 

a gravity gradiometer means; 

a three axis inertial position determining means in- 

cluding a vertical accelerometer means; 

a vertical displacement sensor; and 

an optimal filter means for integrating electronic 
outputs from said gravity gradiometer means, said 
three axis inertial position determining means and 
said vertical displacement sensor. 

2. A navigation system as defined in claim 1 wherein 
said gravity radiometer means measures less than 
every element of the complete local gravity gradient 
tensor and produces estimates of the local gravity field. 

3. A navigation system as defined in claim 1 wherein 
said gravity gradiometer means comprises: 

a plurality of accelerometers mounted on a rotatable 

support mounting means; 

means for rotating said rotatable support mounting 
means; and 

a housing for said means for rotating and said rotat- 
able support mounting means supported by a stabi- 
lized platform. 

4. A navigation system as.defined in claim 3 wherein 

said gravity gradiometer means comprises: 

a shake accelerometer oriented along the axis of said 
rotatable support mounting means and mechani- 
cally coupled thereto; 

a shake actuator mechanically coupled to said shake 
accelerometer; and 

spring bias means for supporting said shake actuator 
and said shake accelerometer within said housing. 

5. A navigation system as defined in claim 3 wherein 
said plurality of accelerometers comprises four acceler- 
ometer. 
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6. A navigation system as defined in claim 3 wherein 
said plurality of accelerometers are mounted tangen- 
tially and equally spaced apart radially on said rotatable 
support mounting means. 

7. A navigation system as defined in claim 1 wherein 
said gravity gradiometer means comprises one to three 
rotating gravity gradiometers. 

8. A navigation system as defined in claim 1 adapted 
to function as a gravity survey system. 

9. A navigation system as defined in claim 1 wherein 
said gravity gradiometer means measures every element 
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of the complete local gravity gradient tensor and pro- 
duces estimates of the local gravity field. 

10. A navigation system, comprising: 

a gradiometer means comprising from one to three 
gradiometers each comprising four rotating accel- 
erometers mounted tangentially at cardinal points 
on a rotating means; EI 

a three axis inertial position determining means in- 
cluding a gravimeter grade vertical accelerometer; 

a depth sensor; and 

an optimal filter means for integrating outputs from 
said gradiometer means; said three axis inertial 


position determining means and said depth sensor. 
* * * * * 


